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Continuous- variable cluster states offer a potentially promising method of implementing a quan- 
tum computer. This paper extends and further refines theoretical foundations and protocols for ex- 
perimental implementation. We give a cluster-state implementation of the cubic phase gate through 
photon detection, which, together with homodyne detection, facilitates universal quantum compu- 
tation. In addition, we characterize the offline squeezed resources required to generate an arbitrary 
graph state through passive linear optics. Most significantly, we prove that there are universal states 
for which the offline squeezing per mode does not increase with the size of the cluster. Simple repre- 
sentations of continuous- variable graph states are introduced to analyze graph state transformations 
under measurement and the existence of universal continuous- variable resource states. 
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I. INTRODUCTION 

Nonstandard models of quantum computation are im- 
portant, both practically and conceptually. On the one 
hand, they lead to new experimental methods to realize 
quantum computers; on the other hand, they offer addi- 
tional insight on the often counterintuitive properties of 
quantum information. Continuous-variable (CV) quan- 
tum computation [T] not only provides a framework for 
description of interacting quantum fields [2] , but also of- 
fers additional realizations of quantum computers when 
each CV mode is assigned a suitable qubit encoding [3J 0]. 
Meanwhile, cluster-state computation [5 showed that the 
implementation of many difficult Hamiltonians may be 
avoided by just applying single-qubit measurements on a 
suitably prepared multi-party entangled resource state; 
hence challenging traditional intuition that the imple- 
mentation of a general unitary operator requires unitary 
evolution. 

CV cluster-state computation is a fusion of these pro- 
tocols [SJ |7] . In addition to its intrinsic conceptual inter- 
est, the formalism presents a potential alternative imple- 
mentation of a quantum computer. Optical CV clus- 
ter states have distinct advantages over discrete ana- 
logues [8J. Any such cluster state may be generated 
deterministically through offline squeezing and passive 
linear optics [3], while all multi-mode Gaussian transfor- 
mations performed through the cluster require only ho- 
modyne detection [TJ. In addition, via alternative tech- 
niques, large CV clusters can be generated in a single 
step using just one optical parametric oscillator (OPO) 
and no interferometer |10j ; some such proposals also have 
significant scaling potential [TTJ [T^J- These features of 
CV cluster states suggest that they offer a fertile experi- 
mental testing ground for the principles of measurement- 
based computation [13]. CV cluster states involving 
four optical modes have been demonstrated experimen- 



tally [HI [33J OS]. 

In this paper, we expand and extend the results given 
in Ref. [7J. First, we apply the CV stabilizer formal- 
ism [TTJ to give simple phase-space and algebraic rep- 
resentations of CV graph states. We then apply these 
tools to compute how graph states transform through 
quadrature measurements and show that there exist uni- 
versal graph states — cluster states — that can be used as 
resource states for the implementation of an arbitrary 
CV circuit 1 . 

Second, we extend the results of Ref. [S] by bounding 
the offline squeezed resources required to construct an ar- 
bitrary graph state to a given precision through passive 
linear optics. These results are applied to several graph 
states of common interest, including linear graph states 
and universal cluster states. We show that the level of 
squeezing required per mode does not grow with the size of 
the cluster state; a necessary criterion to perform quan- 
tum computation efficiently through offline resources. In 
addition, we prove that even if online squeezing is as- 
sumed to be as readily available as its offline counterpart, 
the generation of CV cluster states via offline resources 
remains less costly. 

Third, we detail an explicit optical implementation of 
a non-Gaussian operator through photon counting and 
homodyne measurements and thus propose an explicit 
measurement sequence on CV cluster states that facili- 
tates universal quantum computation. We also present 



In this article, a 'graph state' can have an arbitrary graph, while 
a 'cluster state' must be a member of a family of graph states 
that is universal for quantum computation. The reader should 
be aware that conventions vary in the literature, and these terms 
are sometimes used interchangeably. We will, on occasion, use 
the term 'cluster' on its own, whose meaning at the time should 
be clear from the context. 
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an alternative formalism such that the embedding of non- 
Gaussian resource states allows for universal quantum 
computation entirely by homodyne measurements alone. 
Together these results refine many of the details of the 
CV formalism, offer tools for further development of CV 
cluster-state protocols, and present a variety of poten- 
tially promising and viable experiments. 

The structure of the paper is as follows. Section II 
describes background material on CV quantum compu- 
tation and qubit cluster states that will be useful later 
in the paper. Section III introduces graph states for 
CV modes (qumodes) and describes their stabilizer and 
phase-space representations. Section IV demonstrates 
that such states, coupled with single-qumode measure- 
ments, are capable of implementing any specific unitary. 
Section V explores how CV graph states transform un- 
der measurements and applies these results to construct 
a CV cluster state that may be used as a resource for 
universal quantum computation. In Section VI, the case 
of imperfect CV clusters and the resulting distortions are 
analyzed and discussed. Section VII discusses the optical 
implementation of CV cluster-state computing, including 
the resource requirements for generating arbitrary CV 
clusters and explicit implementation of a nonlinear gate 
that facilitates universal quantum computation. Section 
VIII concludes the paper. 

II. PRELIMINARIES 

In this section, we review some of the background 
knowledge relevant to CV cluster-state computation, 
and its optical implementation. Familiarity with quan- 
tum computation and quantum optics to the level of 
Refs. [HJ [19] is assumed. CV cluster-state computation 
combines the concepts of CV quantum computation and 
cluster states. For a more extensive review of these top- 
ics, please see Refs. [2lfT3]. 

A. Continuous- Variable Quantum Computation 

1. CV State Representations 

In traditional quantum computation, which uses dis- 
crete quantum variables, the basic unit of information is 
the qubit, a system with a two-dimensional Hilbert space 
with computational basis states |0) and |1) and conjugate 
basis states |+) and |— ). The two bases are related by 
the Hadamard operation H. 

The analogue for CV quantum computation is 
the qumode 2 , a quantum system with an infinite- 
dimensional Hilbert space spanned by a continuum of 
orthogonal states \s) q for each s £ K, with orthogonality 



condition (r\ q \s) q — 5(r — s). The conjugate basis states 
are labeled \s) p . The two bases are related by a Fourier 
transform operation: 

I«>p = 4= r dre irs \r) q = F\s) q , 

V J —oq 

1 r°° 

\ s ) q = ^T dre- irs \r) p = Fi\s) p . (1) 

V J — oo 

The unitary operator F is defined by this relation. In 
quantum protocols, qumodes may be used to encode 
qubits (e.g., the GKP encoding [4 or a coherent-state 
encoding [5]), or they may be employed directly for CV 
quantum computation [TJ [20] . 

We may now define corresponding observables, posi- 
tion q and momentum p, such that q\s) q — s\s) q and 

P \ s )p ~ s l s )j)' with [q,P] — i where % = 1. Here, p is the 
generator of positive translations in position, while — q 
is the generator of positive translations in momentum. 
Thus, we can write an arbitrary position and momentum 
eigenstate as 

\s) q = X(s)\0) q , \s) p = Z(s)\Q) p , (2) 

where X(s) — e~ ls ^ and Z(s) — e ts ^ represent displace- 
ments in the computational and conjugate basis, respec- 
tively. An arbitrary pure quantum state \<p) of a CV 
system may be decomposed as a superposition of either 

\ s ) p or \ s )„- 

While the computational basis or its conjugate is un- 
countable, any physical state \(f>) may nevertheless be de- 
composed into a count ably infinite basis. For particles 
in a harmonic trap or quantum optical fields we can use 
the Fock basis of definite particle number {|0) , |1) , . . .} 
where h = a) a is the number operator, with n\n) — 
n\n), the usual bosonic commutator [a, a^] = 1, and 
a = (q + ip) j In the terminology of quantum optics, 
q and p are referred to as the 'position quadrature' and 
'momentum quadrature' for a given mode, respectively. 

A qumode is in a minimum uncertainty state if the 
product of the quadrature deviations Aq and Ap is min- 
imized, i.e., AqAp = The ground or vacuum state 
|0) defined by a |0) = is an example of particular the- 
oretical and practical interest and represents a Gaussian 
superposition centered about in either the computa- 
tional or the conjugate basis: 

10) = / ds ^ ^ = / dse ~ s2/2 ^r> ■ 

The vacuum state is a specific example of a Gaussian 
state whose quadratures exhibit Gaussian statistics. 

The state of a single qumode can be described by its 
Wigner Function [5T] 

W(x,y) = ± Jdw(x~^\^\x+^)^y. (4) 



2 We use the terms 'mode' and 'qumode' interchangeably. 



The Wigner function is a useful tool for describing arbi- 
trary Gaussian states, which are completely determined 
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by the first and second moments of the quadratures. Any 
state with a Gaussian Wigner function is, by definition, a 
Gaussian state. For instance, the Wigner function of the 
vacuum state is e~( x +y V 71 "; a multivariate Gaussian dis- 
tribution with a variance of 1/2 in both quadratures. A 
multi-mode state such as a CV cluster state is described 
by a multi-mode Wigner function, a straightforward ex- 
tension of Eq. Q . Multi-mode Gaussian states are then 
given by a second-moment covariance matrix and a first- 
moment vector [2J. 

2. Gaussian Transformations 

In quantum optics, the Hamiltonians corresponding 
to the experimentally most feasible interactions are at 
most quadratic in q and p. Such interactions trans- 
form Gaussian states to Gaussian states, and are re- 
ferred to as Gaussian transformations or linear unitary 
Bogoliubov (LUBO) transformations. If we collect the 
quadrature operators into an operator-valued vector v = 
(<?i> <Z2) • • • >Pi)P2i ■ ■ -) T j then a general Gaussian transfor- 
mation U transforms v according to: 

U ] vU = Lv + c, Det(i) = l, (5) 

where L is a In x 2n symplectic matrix and c is a vector 
of 2n constants that represent quadrature displacements. 
We list a number of standard single-mode Gaussian trans- 
formations that will be used in this paper, along with 
their associated Heisenberg action on the quadrature op- 
erators. 

(a) Rotations: R{6) = e' e (3 2 +p 2 )/2 rotates a state coun- 

terclockwise in phase space by an angle 9. Rota- 
tions are also referred to as phase shifts. 

(iH™rr/)(JH*(<->(!) «> 

where Mr{6) is the rotation matrix describing the 
linear Heisenberg action on the quadrature opera- 
tors. Note that R(n/2) = F. 

(b) Quadrature Displacements: Z(s) — e ts ^ dis- 

places a state in phase space by s in momentum. 

(IMfM°) 

Similarly, X(s) = e _lsj3 displaces a state in phase 
space by s in position. Note the sign in the expo- 
nential of each. 

(c) Squeezing: S(s) = e -*M»)(li»+fi$)/2 sque ezes the 

position quadrature by a factor of s, while stretch- 
ing the conjugate quadrature by 1/s. 




where Ms(s) is the squeeze matrix describing the 
linear Heisenberg action on the quadrature opera- 
tors. 

(d) Shearing: -D2,<j(s) = e IS * / 2 shears a state with re- 
spect to the q axis by a gradient of s. The shearing 
operator e lsil I 2 is also referred to as the phase gate. 

({)-(!!)(«)-**>(*) 

where Mn(s) is the shearing matrix describing the 
linear Heisenberg action on the quadrature opera- 
tors. 

Operations (a) and (b) correspond to the most readily 
available single-mode Gaussian transformations, requir- 
ing only phase shifts and coherent-state sources. To ac- 
cess all possible single-mode Gaussian transformations, 
we will need squeezing interactions to stretch and com- 
press phase space uncertainties. Two such operations 
are given by (c) and (d). In experimental quantum 
optics, such interactions require nonlinear optical pro- 
cesses (while the Heisenberg in-out relations remain lin- 
ear). Typical methods involve optical parametric am- 
plification, which allows one to generate squeezed vac- 
uum states S(s) |0). We refer to such processes as offline 
squeezing, solely involving the preparation of squeezed 
vacuum states. 

Offline squeezing contrasts with the online squeezing, 
where the squeezing operator is applied "on-line" to an 
arbitrary state of the electromagnetic field. In experi- 
mental quantum optics, it is common to refer to S(s) |0) 
as a state with 101og(s 2 ) dB of squeezing, alluding to 
the view that squeezing can be regarded as a physi- 
cal resource [2J. While the generation of reasonably 
high levels (10 dB) of offline squeezing can be experi- 
mentally achieved |22j . online squeezing |23j is far more 
demanding, and is currently only experimentally viable 
for modest values of s, for instance, online squeezing of 
2.5 dB [23] utilizing offline squeezed ancilla states [2"S"] . 

An arbitrary single-mode Gaussian transformation 
may be decomposed into (a) rotations, (b) quadrature 
displacements, and either (c) squeezing or (d) shear- 
ing operations. The addition of a two-mode Gaussian 
gate, such as a beamsplitter or Cz = e 1 ^®^, allows for 
arbitrary multi-mode Gaussian transformations. To ac- 
count for imperfect Gaussian transformations, e.g. af- 
fected by photon losses and thermal excess noises, Gaus- 
sian unitary gates are generalized to Gaussian operations 
(Gaussian completely positive maps) |26j . These also in- 
clude Gaussian measurements such as homodyne detec- 
tion. Any quantum evolution consisting solely of Gaus- 
sian operations on Gaussian states may be efficiently sim- 
ulated on a classical computer [20] . Therefore, some sort 
of non-Gaussian element is required for universal quan- 
tum computation. In fact, at least in principle, any such 
element will do [1 . 
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3. A Universal Gate Set 

We follow the definition of universal CV quantum com- 
putation outlined in Ref. [Tj . A system is universal if it 
can simulate the action of a Hamiltonian consisting of a 
general polynomial of p and q to any fixed accuracy. 

For a single qumode, all Gaussian operations together 
with any single nonlinear (non-Gaussian, at least cubic) 
interaction are capable of universality [Ij. For example, 
the set of gates Dk, q (s) — exp(isq k /k), for k = 1,2,3 
for all s E R, together with the Fourier transform F, 
are sufficient for universal single-mode quantum compu- 
tation (that is, this set can be used to implement any 
single-mode unitary operation up to arbitrary accuracy) . 
Here -Di j( j(s) is a displacement, D2,q(s) is a shear, and 
D 3 ^(s) is the cubic phase gate [I]. Adding to this set 
any nontrivial two-mode interaction allows for universal 
quantum computation. For theoretical simplicity, here 
we shall use the Cz gate (defined above) to complete the 
universal set, while another possibility is a simple beam- 
splitter interaction. 

It should be noted that such statements about uni- 
versality do not account for noise. Presently, all gen- 
eral error correction codes require discretization at some 
level. Hence, currently CV quantum computation is only 
known to be possible for discretized encodings of CVs. 

B. Cluster-State Computation 

In the qubit-based cluster-state model jS], quantum 
computation is implemented by a series of single-qubit 
measurements on a specially prepared, entangled state 
of many qubits, most generally, referred to as a graph 
state [27]. Such states may be conveniently described 
by graphs. A graph G = (V, E) consists of a vertex set 
V = {«i}™ =1 and a set of edges E. We say that two 
vertices, vt and Vj, are neighbors if there exists an edge 
(vi,Vj) € E that connects them. For an introduction to 
graphs, see Ref. [28] . 

For any undirected, unweighted graph G = (V,E) hav- 
ing no self-loops, we can construct a corresponding graph 
state as follows. For each vertex of G, we initialize a qubit 
in the state |+) = ^(|0) + |1)). For every edge in G link- 
ing two vertices, we apply a CSIGN gate (which is some- 
times called the CPHASE gate) to the two corresponding 
qubits. Any unitary operation can be implemented on a 
tailor-made graph state using an appropriate sequence of 
single-qubit measurements. 

The stabilizer formalism |29j offers an efficient way to 
represent any graph state. A state \<f>) is stabilized by an 
operator K if it is an eigenstate of K with unit eigenvalue, 
i.e., K \4>) = \<f>). The set of stabilizers form an abelian 
group under operator multiplication. If such a set exists 
for a given state, then we call that state a stabilizer state, 
and we may use the generators of its stabilizer group to 
uniquely specify it. The stabilizers for qubit graph states 
are well known. Given that \<fi) is an n-qubit graph state 



with associated graph G = (V,E), it is stabilized by 

Ki = Xi [] z i (8) 

jeN(i) 

where N(i) denotes the set of indices that define the set of 
vertices that neighbor v i} i.e., N(i) = {j | (vj,Vi) € E}. 
The operators X and Z are the usual Pauli operators for 
qubits. 

There exist universal families of graph states that 
may be used to implement any unitary operation solely 
through the choice of single-qubit measurements made 
on it. Originally, such states are called cluster states, 
and cluster-state computation involves the implementa- 
tion of arbitrary algorithms solely by an adaptive local 
measurement scheme. The scheme involves only single- 
qubit measurements and is called 'adaptive', because the 
choice of measurement bases depends both on the algo- 
rithm to be implemented and, in general, on the mea- 
surement outcomes during a cluster computation. Clus- 
ter states, when combined with adaptive local measure- 
ments, are thus universal resources for quantum com- 
putation [5]. For more recent developments on possible 
resource states for universal quantum computation and 
their requirements, see Refs. [30 ] l3T ] 132 ] [33] . 



III. CONTINUOUS- VARIABLE GRAPH STATES 

The concepts of qubit cluster-state computation can be 
extended to the continuous-variable regime. We outline 
CV graph states [B], which can be used as resource states 
for universal CV quantum computation [TJ. We then in- 
troduce milliners, a variation of the CV stabilizer for- 
malism [S1[T7], and use them to compute how CV graph 
states transform under quadrature measurements. 

The basic premise of CV graph states may be ob- 
tained by replacing elements of qubit cluster-state com- 
putation with their CV analogues: |+) becomes |0) , X- 
measurements are replaced by measurements of p (and 
Z with q), and the CSIGN interaction is replaced by the 
Cz = e l ^^ j gate, which is used to entangle nodes i and j. 
Each CV graph state can also be defined by a graph 
G = (V,E), where the set of vertices V corresponds to 
the individual qumodes, and the edge set E determines 
which qumodes interact via the Cz operation. 

It should be mentioned that one way to generalize the 
idea of a CV graph state is to use weighted edges for 
the graph. The edge weights specify the strength of the 
Cz interaction between the connected nodes: Cz(t) = 
e »*9®9 j where t is the edge weight. CV weighted graph 
states have a variety of uses jlOj [34j [35] , but in this article 
we will confine further discussion to unweighted graphs 
(or, equivalently, graphs whose edge weights are all +1). 



A. Stabilizers and Nulliflers 

Analogous to the case for qubit graph states, the sta- 
bilizer formalism for CV systems [HI [T7] can be used to 
specify any CV graph state completely [M] . We say that 
a zero-momentum eigenstate |0) is stabilized by X(s) 
for all s, since it is a +l-cigenstate of those operators. 
This holds even though X(s), being non-Hermitian, is 
not an observable. Notice that if K stabilizes \4>), then 
UKW stabilizes U\(f>). This observation, together with 
the relation e lqiq2 pie~ lqiq2 = pi — $2, allow us to write 
the stabilizers for an arbitrary CV graph state \<fi) on n 
qumodes with graph G = (V,E): 



(a) 



K i (s)=X i (s) J] Zjis), 
jeN(i) 



1, 



(9) 



for all s G I, where N(i) is defined as before in Eq. ([8]), 
and the subscript indicates which qumode the displace- 
ment acts on. 

This group is conveniently defined by its Lie algebra, 
the space of operators H such that H \ <j>) = 0. We refer 
to any element of this algebra as a nullifier of \<p) and 
the entire algebra as the nullifier space of \4>). Being 
Hermitian, every nullifier is an observable. Any ideal 
graph state has a simple nullifier representation. 

Theorem 1 The nullifier space of an n-qumode graph 
state \4>) with graph G — (V, E) is an n- dimensional vec- 
tor space spanned by the following Hermitian operators: 



Hi =p i- ^ q 
jeN(i) 



i = 1, 



(10) 



That is, any linear superposition H = CiHi satisfies 
H \<j>) = 0. Note that [H h Hj] = for all (i,j). 

Proof: Every stabilizer from Eq. ([9| is the exponential 
of a nullifier in this space. Specifically, Ki(s) = e~ lsHi 
for all s G M., with i = 1, . . . , n. M 

In Fig. [l] we illustrate this formalism. Note that the 
nullifier space for a given state does not have a unique set 
of milliners (a basis) that defines it, since linear combi- 
nations of milliners give another nullifier. Nevertheless, 
Eq. (10) is a standard set that can easily be obtained 



from a given graph. 



B. Wigner Representation 

The Wigner function can be useful as an extension 
to the nullifier formalism. It encapsulates the simplic- 
ity of the nullifier formalism, while maintaining the intu- 
ition afforded by an explicit representation of the state, 
and importantly, continues to be useful for non-ideal CV 
cluster states. Since the arguments of a Wigner func- 
tion behave identically to the milliners under Gaussian 
transformations, they may also be written by inspection. 





(i+3j+3) 



0+3j) 



FIG. 1: Nullifiers give an efficient description of ideal graph 
states. The nullifier space of the linear graph state on three 
nodes (a) is spanned by pi — 52, P2 — qi — <73 and ps — qa. 
The infinite two-dimensional lattice (b) is nullified by Hij = 
Pi j —qi-i,j —Qi+l,j for each coordinate (i, j). 



The Wigner function of an ideal n-qumode graph state 
\4>) with graph G — (V,E) is a function of 2n variables 
on the scalar- valued vectors q = (q±, . . . , q n ) and p = 
{Pi, ■ ■ ■ >Pn)- Explicitly, 



W(q,p) = l[e( qi )S(H t ) 



(11) 



where Hi, i — 1, ... ,n are the standard nullifiers of \<p) 
(with each of the operators cji and Pi replaced by scalar 
variables qi and pi, respectively), 6(x) is the Dirac delta- 
distribution, and e(x) is the infinite uniform distribu- 
tion. Ideal CV graph states are highly singular, so for 
all practical purposes, S(x) and e(x) should be consid- 
ered limits of a normalized Gaussian whose variance, re- 
spectively, vanishes and extends to infinity. For example, 
the Wigner function of Fig. [TJa) is e(9i)e('?2)e((J , 3)^(pi — 
q 2 )S(p2 - qi- q 3 )S(p 3 - q 2 ). 

Wigner functions can also be used to define an ex- 
tended class of generalized graph states. Whereas an 
ideal graph state with associated graph G may be de- 
fined by the action of appropriate Cz gates on n mo- 
mentum eigenstates, a generalized graph state replaces 
each momentum eigenstate with some arbitrary quantum 
state \4>i). If \4>i) has a corresponding Wigner function 
Wi(qi,pi), then the Wigner representation of the result- 
ing generalized graph state is given by 



W(c l ,p) = l[W i (q i ,H i 



(12) 



j=i 



Such states are used extensively when we perform compu- 
tations with graph states and when we extend the graph 
state formalism to realistic situations where momentum 
eigenstates need to be approximated. 
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IV. QUANTUM COMPUTATION ON CV 
GRAPH STATES 



CV graph states are a resource for CV quantum com- 
putation. For any given CV unitary U, and any given 
input \<p), there exists an appropriate graph state such 
that by entangling the graph state locally with \<j>) and 
applying an appropriate sequence of single-qumode mea- 
surements, U\4>) is computed. 

To justify this statement, we first show that there 
exists a |</>)-dependent quantum state on a system of 
qumodes that collapses into U\<f>) (modulo known sin- 
gle qumode Gaussian operations) when an appropriate 
sequence of single-qumode measurements is applied. We 
then demonstrate that this ^-dependent quantum state 
can be efficiently constructed using an appropriate graph 
state as a resource. 



A. Measurement-Based CV Quantum 
Computation 

To implement any unitary operation on k qumodes, we 
apply the following algorithm. We first introduce a graph 
G = (V,E). We designate k vertices of G as input ver- 
tices, and another (possibly overlapping) set as output. 
Call these sets V\ n and V^ ut . The following algorithm 
computes U \ <f>): 

1. The qumodes corresponding to the vertices in V{ n 
encode the input state \4>), while the qumodes cor- 
responding to the other vertices are each initialized 
in the state \0) p . 

2. For each edge (vj,Vk) € E, apply Cz = e^^ k be- 
tween vertices j and k. Since all Cz operations 
commute, their order does not matter. 



3. Measure each vertex m, for all Vi (jL V wt in a basis 
of the form M, = e - ift ^pe ifi ^ , where /;(£) is, 
in general, a polynomial of q. The exact form of 
each fi is dictated by the unitary we wish to im- 
plement and the result of measurements on prior 
modes. Without loss of generality, we can label the 
vertices such that they are measured in numerical 
order. 

4. The remaining unmeasured qumodes encode U\<j>), 
modulo known single-mode rotations and transla- 
tions. 

The above algorithm may be implemented by using an 
appropriate graph state as a resource. This algorithm 
is universal. Given any unitary U, there always ex- 
ists an appropriate graph G = (V, E) and designations 
Vi n , V out C V such that the above algorithm implements 
U. 



B. Proof of Universality 

To prove the above procedure is universal, we need 
to show that it can implement (a) single-mode Gaussian 
operations, (b) the cubic phase gate, and (c) the Cz gate. 
First observe that (c) may be implemented trivially by a 
two vertex graph where both vertices are designated as 
both input and output. No measurements are required. 

The implementation of (a) and (b) also each involve a 
two- vertex graph. We designate one vertex as input and 
the other as output. Consider first the case where the 
input mode is measured in the p basis: 



p j= m 
X{m)F\4>) 

Given an input \(f>) ® |0) = J dsf(s) \s) |0) , the state 
of the system after the application of the Cz gate is 



C. 



:(W®|0> P )= J dsf{s)\s) q \s) p . (13) 



Measurement of p on the first mode with associated result 
to, as shown, collapses this state to 



/ out 



oc J dsf{s)({m\ p \s)^ \s) p 
ex j d S f{ S )e~* sm \s) pl 

dsf{s)\s) p = X{m)F\<t>). (14) 



The effect of this circuit is to apply the identity opera- 
tion, modulo a known quadrature rotation and displace- 
ment. Obviously, a transformed input state of D q \<f>), 
for any Dq = e l $^ diagonal in the computational ba- 
sis, would result in output X(m)FD q \<j>). However, this 
same transformation can be induced by an appropriate 
measurement. We can see this immediately by writing 
out the associated circuit. Since D q commutes with Cz, 
the circuit 



10), 



D S 



)= TO 



X{m)FDq 



must have the desired output. The application of Dq 
followed by a p-measurcmcnt has an identical effect to a 
measurement in the rotated basis Pf( q ) = DtpDq. Hence, 
any unitary diagonal in the computational basis can be 
implemented by a single measurement of Pf(q)- Measure- 
ments on two qumodes of a three-qumode cluster equates 
to a repeated application of this circuit, resulting in the 
output 



X{m 2 )FDqX{ mi )FD q 



= X(m 2 )FX( mi )D (q 



(q+rn^FDq \<j)) 

X{ m2 )FX{ mi )FD ( _ P+mi) D q \4>) i (15) 
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where D$ = e l ^ v ^ for any operator v, and we have used 
the relations 



X(—m)qX(m) — q + m, 
Z \—rn)pZ (m) = p + m. 
F\-q)F=p, 
F^pF = q, 



the last two of which give 



F ] Z{m)F = X(m), 
F i X(m)F= Z{-m). 



(16) 
(17) 
(18) 
(19) 

(20) 
(21) 



If, instead of Pfif\ as shown above, we had mea- 
sured the second mode in the outcome-dependent ba- 
sis fiff_q_ mi \, the output would instead be |0) out = 
X(m2)FX(mi)FDpDq \<fi). Thus, the ability to measure 
the second mode in the basis pj(_^_ mi ) allows determin- 
istic implementation of Dp. 

By concatenating this circuit a sufficient number of 
times, any single-mode operation can be implemented de- 
tcrministically by alternating application of Dq and Dp 
(with generally different Ds each time) [I]. Note that the 
measurements required to implement these gates (beyond 
the first one) are necessarily adaptive — that is, our choice 
of the measurement basis is generally dependent on previ- 
ous measurement results. Also notice that the final result 
is modified by a measurement-dependent Gaussian oper- 
ation (X(m2)FX(mi)F , in the simple case illustrated). 
This need not be corrected. As long as we keep track of 
it, it can instead be considered as just a change of basis 
for the final state. 

Another useful way of approaching the question of uni- 
versality in the CV context is to consider implementing 
Gaussian operations and then, separately, non-Gaussian 
operations. Single-mode Gaussian operations require the 
ability to implement e lsq (quadrature displacements) and 



(shears) for all set, plus the Fourier transform F. 



J.sq/2 

The Fourier transform is obtained for free simply 
through the Gaussian correction applied with each mea- 
surement. Displacements e lsq are easily implemented, as 
well: just measure p sq = p + s, which is the same as 
measuring p and adding s to the result 3 . In this case, 
dependence on previous measurement outcomes is trivial 
since p s (q+m) is also equal to p+ s, and thus, no adapta- 
tion is required at all. 

Shearing transformations e %s ^l 2 correspond to mea- 
surements in the basis p S qi /2 = P + s <7- In the 
case that adaptation for previous measurements is re- 
quired, the new measurement basis would be of the 
form p s (q+ m ) 2 ii — p+sq+ms, which differs from the orig- 
inal basis only by a measurement-dependent constant. 
This can be accounted for trivially by measuring in the 
original basis and adding ms to the result. 



Since the "adaptation" required for previous measure- 
ment outcomes is trivial for all measurements necessary 
to implement Gaussian operations, such measurements 
may be made in any order — or simultaneously. This 
property is known as parallelism [7\. Qubit cluster-state 
computation has an analogous property with the same 
name, whereby measurements that implement Clifford 
group operations can be performed in any order |13j . 

The above measurements allow for any Gaussian op- 
eration to be implemented. But universality requires 
non-Gaussian operations as well [20 . The cubic phase 
gate e lsq / 3 allows implementation of all single-mode non- 
Gaussian operations [Tj and may be implemented via a 
measurement in the basis p sq -3 / 3 = p + sq 2 . The dif- 
ference with this gate is that when an adaptive imple- 
mentation is required, the physical basis is now differ- 
ent: p s (g+ m )3/3 = p + sq 2 + |msq + |m 2 s, due to the 
presence of the noncommuting m-dependent term ^msq. 
Accounting for this difference requires physically chang- 
ing the basis of measurement (unlike the last term ^m 2 s, 
which can be eliminated simply by shifting the result). 
As with qubit quantum computation, a general CV quan- 
tum computation will require adaptive measurements for 
the non-Gaussian (non-Clifford) part of the computation. 
What these measurements are in an experimental context 
will depend on the choice of the physical implementation. 
In section VII, we propose one possible method that uses 
photon counting. 

A sequence of singlc-qumode unitaries and wires is im- 
plemented by a sequence of measurements on a linear 
graph state. Cz gates are implemented by edges be- 
tween linear clusters (Fig. [2]). Thus, we may apply the 
algorithm described to implement any given CV unitary 
on an arbitrary input state. This proves universality. 



C. Graph States as Resources 

Observe that steps 1 and 2 generate a special class of 
generalized CV graph states. Each of the input qumodes 
arc initially set to encode the inputs of the desired quan- 
tum computation, rather than the standard momentum 
eigenstates. The resulting cluster has the Wigner repre- 
sentation 

W(q,p) =l[W i (q i ,H i )l[e(q j )6(H j ), (22) 
iei j&i 

where X = {i \ Vi S Vin} is the set of indices that corre- 
sponds to the input qumodes. 

To see that the standard graph state with graph 
G(V, E) may be used as a resource for the algorithm, 
we show that it may be used to efficiently generate 
clusters of the form given by Eq. (22 1. Let the input 



There is a sign error in the corresponding derivation in Ref. 0. 



state be initially encoded on k qumodes, which we label 
{til, U2, . . . , Uk}, and the input vertices of the graph state 
given by V\ n = {v±, V2,---, We proceed as follows: 

1. Apply a Cz operation between each qumode pair, 
(ui, Vi), for each i = 1, . . . k 
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To prove the existence of such universal resources, we 
explore how CV graph states transform under single- 
mode quadrature measurements. These tools are then 
applied to show that there exists a universal CV graph 
state, which, when appropriate quadrature measure- 
ments are applied, collapses to the specific CV brickwork 
state that implements any given quantum circuit. Such 
universal graph states are called CV cluster states. 



(b) 




A. Graph State Transformations 

The nullifier formalism is ideal for computing how 
graph states transform through quadrature measure- 
ments. In this formalism, we describe a measurement pi 
on the i th qumode, with measurement result mi, as fol- 
lows. We first rewrite the milliners in a basis such that 
only one element, say Hi, does not commute with our 
basis of measurement. Hi is then replaced with pi — m,, 
and in all other milliners, pi is replaced with m*. Mea- 
surements in the fa basis are treated analogously. The 
details of this formalism are outlined in Appendix |A| 



FIG. 2: (a) Any unitary on multiple qumodes may be written 
as a quantum circuit consisting of Cz and single-qumode uni- 
taries diagonal in either the position or momentum basis [T]. 
(b) Any such circuit may be directly translated into an appro- 
priate graph state. Here the arrowed qumodes are measured 
in the appropriate basis that implements their corresponding 
single-qumode unitary. All other non-output qumodes are 
measured in the p basis. 



2. Measure each of the modes Ui, resulting in mea- 
surement results rrii, i = 1, . . . , k. 



The resulting cluster state is identical to Eq. (22 1, modulo 



known single-mode quadrature displacements and rota- 
tions on each Vi that can be accounted for throughout 
the remainder of the computation. 

Thus, the circuit of Fig. [2ja) may be implemented by 
using a standard graph state with the graph shown in 
Fig.[2^b). We refer to such graph states as CV brickwork 
states, alluding to similar results in qubit cluster-state 
computation |36j . 



V. UNIVERSAL CLUSTER STATES 

So far, we have discussed the construction of specific 
graph states tailored for a specific quantum computation. 
Like qubit clusters, there exist classes of universal CV 
graph states that may be used to implement an arbitrary 
CV operation. This is of more than theoretical interest, 
since it facilitates the development of physical systems 
that are tailored to generate one particular state. This 
state can then be used as a universal resource. 



1. Vertex Removal 

A computational-basis measurement on a qumode re- 
moves it, along with all edges that connect it, from the 
cluster. Consider a measurement fa on the i th mode of 



\<f). Equation (10 I indicates that Hi is the only noncom- 
muting basis element. Therefore a measurement with re- 
sult rrii transforms Hi into fa—mi and replaces fa with Wj 
in all others. Explicitly, for each j =/= i, 



H 3 H 3 „ 

qi — >rrii 

(23) 

The resulting state corresponds to the graph state of G 
with vertex Vi removed, modulo known quadrature dis- 
placements. This operation is useful for creating a CV 
graph state that corresponds to the subgraph of some 
original resource state. In addition, it functions as a 
handy "delete" button, and can be used to "amputate" 
corrupted parts of a cluster state. To summarize, if \<p) 
is the graph state defined by a graph G = (V, E), a q- 
measurement on a mode i results in the graph state with 
associated graph G\{vi}, modulo known corrections, i.e., 
the graph resulting from removal of vertex Vi along with 
all edges connecting to u,. Thus, a computational mea- 
surement removes a given node from the cluster. 
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FIG. 3: Any CV graph state may be generated by appro- 
priate single-mode measurements. Computational-basis mea- 
surements (blue orbs) remove unwanted nodes. Momentum- 
basis measurements (green orbs) are then employed to shorten 
the "wires" within the cluster. 



appropriate graph state for simulating any given circuit 
(Fig. [3j. The graph state that corresponds to a planar 
square lattice is thus a resource for universal CV quan- 
tum computation and is therefore a CV cluster state. 

In practice, of course, lattices are always of finite size, 
just as are all quantum circuits. Therefore the complex- 
ity of the quantum computation one wishes to perform 
is constrained by the size of the original resource state. 
Since the size of the required cluster grows linearly with 
the number of fundamental one- and two-qumode gates 
and also grows linearly with the number of qumodes, CV 
cluster-state computation is efficient. As in the case of 
qubits, any algorithm of polynomial gate complexity can 
also be implemented by a resource state of polynomial 
size. 



VI. THE EFFECTS OF FINITE SQUEEZING 



2. Wire Shortening 

Sometimes we may also wish to remove a given vertex 
but still preserve the connections of its neighbors. This 
transformation is useful, for example, when we wish to 
shorten linear graph states. Consider p-measurements on 
the two inner nodes of the four-node linear graph, which 
has nullifier basis 

{Pi - fa,P2 -qi- fa, fa - fa- fa,P4 - fa}- (24) 

Since we will be making measurements of p2 and P3, we 
want a new basis in which only one nullifier fails to com- 
mute with f>2 and only one other fails to commute with pV 
We construct this basis out of linear combinations of the 
original basis elements, resulting in 

{Pi - fa,P2 -P4- fa,-Pi+P3 - fa,P4 - fa} (25) 

Measurements of p2 and P3, with outcomes m-i and 7713, 
respectively, collapse the cluster into a new graph state 
with milliners {7712 — p4 — fa , 7713 — p\ — fa}. This is equiv- 
alent to the graph state of the two-qumode cluster, mod- 
ulo a known quadrature displacement and a reflection in 
phase space about one of the nodes. Thus, measurements 
in the momentum basis allow us to effectively "shorten" 
linear graph states. 



B. The Universal Resource State 

Recall that any quantum circuit may be implemented 
by a sequence of measurements on a specifically tai- 
lored brickwork state (See Fig. [2jb) ) . The graph for 
such states is always a subgraph of a sufficiently large 
two-dimensional square lattice (see Fig.[ljb)). The two 
transformations outlined above allow us to carve out an 



The ideal framework of CV quantum computation in- 
volves the use of momentum eigenstates. Such states 
cannot be normalized and are thus an idealized abstrac- 
tion. Any practical implementation must necessarily ap- 
proximate these states. One way to do so is by replac- 
ing each zero-momentum eigenstate with a vacuum state 
that has been finitely squeezed in the momentum quadra- 
ture. In this section, we detail the resulting distortions 
imposed on any quantum computation that uses cluster 
states made from these approximate states. 

Suppose we use states of finite squeezing, i.e., S{s) |0) 
(where |0) represents the vacuum) for some large s, in 
place of momentum-quadrature eigenstates. The result- 
ing graph state obtained will not be ideal. Formally, 
we say that the resulting graph state \4>{s)) is of accu- 
racy s. Such states are generalized CV cluster states, 
and Eq. (12) allows us to write down their Wigner rep- 



resentation: 



W(q,p) = l[G s (q j )G 1/s (H j ), 



(26) 



where G s (q) = (7rs 2 ) -1 / 2 exp(— q 2 /s 2 ) represents a Gaus- 
sian distribution with variance s 2 /2, and {Hj} are the 
milliners of \(f>(s)) from Eq. (10 1. Thus, these general- 
ized CV cluster states are Gaussian states. Observe that 
G s converges to a uniform distribution and G\/ s to a 
(5-peaked distribution in the limit of large s, in agree- 
ment with the Wigner function for an ideal graph state, 
Eq. O. 



To analyze the resulting distortions, we first consider 
the special case of simple state teleportation, where only 
p-measurements are made and the input state propagates 
through the cluster without any intended manipulation. 
This result may then be extended to arbitrary measure- 
ments and the implementation of universal gates. 
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A. Distortions in State Propagation 

Consider the state resulting from a p-mcasurcmcnt on 
a 2-qumode cluster, with input state \<f>) specified by the 
Wigner function W m {q,p). In terms of the circuit model, 
this is represented by 



_pj= m (27) 

S(s)\Q) -I X(m)F\4,)' 

)' is a distorted version of \d>), which will be 



where 

specified below. We can analyze this circuit as follows. 
After entangling the input, the state of the system is 
given by 



W in (qi,pi - g 2 )G s (g 2 )G 1/s (p2 - qi) 



(28) 



A pi-measurement with outcome m yields the output 
state 



P(m)W out (q,p) = G s (q) J drW- m {T,m - q)G 1/s {p - r), 
= G s {q) [(W ln * 1 G 1/s )(p,m-q)] , 

(29) 

where *i denotes a convolution with respect to the first 
argument of W, and P{m) is the probability of measure- 
ment outcome m. P{m) multiplies the resulting, normal- 
ized pure state W ut to give the actual expression on the 
right-hand side. 

What we would like to know from this toy example 
is how the imperfect squeezing affects the encoded state 
under the cluster-state implementation of the identity 
gate. A good way to see this effect is to undo the unitary 
correction X(m)F and compare the result W{ n to the 
original input state W- ln : 

P(m)W{ n = G a (m-p) *i G 1/s )(q,p)] . (30) 

The Wigner function W( n (q,p) corresponds to \<f>) in Cir- 
cuit (27 1. This means that with respect to the quantum 



information to be teleported, the Gaussian envelope is 
dependent on the measurement outcome to. Some values 
of m will result in an envelope that overlaps the (nonneg- 
ligible) support of W- m , while other more extreme values 
of to will result in a strongly shifted envelope that cuts 
off large portions of the support of W[ n . Thus, the actual 
success of any instance of teleportation depends strongly 
on the measurement outcome to. 

On the other hand, we can instead talk about the aver- 
age state (a mixed state) that results from teleportation 
when we average over all possible measurement results m 
using their corresponding probabilities P(m). This state 
is easily calculated using Eq. (30 1: 



W avg = {W{ n ) = / dmP{m)W( n = W in *i G 1/s . (31) 

Thus, the average effect on the quantum information due 
to teleportation using finitely-squeezed resources is just 



the addition of a variance of l/(2s 2 ) noise units on the 
q-quadrature. Repeated application gives us the result- 
ing average distortion when a chain of p- measurements is 
used to teleport an initial state W- ln down a linear cluster: 



avg — Wi n *1 Gi/ S * 2 Gi/ S *i Gi/ 



s *2 • ' 



(32) 



In summary, when propagating quantum information 
through a chain of finite accuracy s, in every single shot, 
pure, conditional output states are created with Gaus- 
sian envelopes applied to the input state in alternating 
quadratures and with the measurement results {to^} de- 
termining their respective centers. More typical, when 
CV quantum information is teleported through a chain 
of finite accuracy s, on average, l/(2s 2 ) units of noise are 
added alternately between the two quadratures resulting, 
in general, in a mixed output state. Whether a single- 
shot or an average picture is applicable depends on the 
actual experimental implementation and the encoding of 
the signal states. 



B. Distortions in Universal Gate Teleportation 

The distortions derived above, caused by finite squeez- 
ing, apply to all singlc-qumode measurements. To sec 
this, consider the application of an arbitrary singlc- 
qumode unitary D, diagonal in the computational basis, 
by measuring in the D^pD basis. Since D and Cz com- 
mute, this is equivalent to standard teleportation with 
input D i.e., 



D\4>) 
S{s) |0> 



X{m)F{D\4>))' 



Thus, the resulting output state is again the expected 
output state in the limit of ideal graph states, subjected 
to the distortion given by Eq. (30 1. Therefore, the use 



of finite squeezing results universally in the addition of 
Gaussian noise that "blurs out" the details in the momen- 
tum and position quadratures, alternating between them 
at each step. The magnitude of this noise is inversely 
proportional to the accuracy of the cluster and grows 
linearly with the length of the cluster. This noise can 
potentially be reduced by the use of redundant rails [9 . 
However, such redundant, multiple-rail encoding requires 
a larger amount of squeezing resources for creating the 
corresponding graph state. We will get back to this point 
in the following sections on optical cluster-state genera- 
tion and computation. 



VII. OPTICAL CLUSTER STATES 

The optical implementation of CV cluster states holds 
particular promise and features a number of advantages 
over its discrete- variable counterpart [8j [37] . With opti- 
cal qubit cluster states, the entangling operation that is 
used to generate a cluster state is highly nonlinear, and 
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its proposed implementations are nondeterministic. This 
results in significant overhead and presents an impedi- 
ment to the generation of large-scale clusters. 

While challenging nonlinear operations are still re- 
quired for universal quantum computation, the genera- 
tion of CV cluster states with current technology is en- 
tirely deterministic. In particular, 

1 . Any CV graph state can be prepared completely via 
the interaction of squeezed vacuum states through 
a network of linear optical elements. Not only do 
we avoid the need for nonlinear interactions, but 
online squeezing is also unnecessary. 

2. Once the cluster state is prepared, any multi- 
qumode Gaussian operation may be implemented 
entirely by quadrature measurements (homodyne 
detection) . 

3. The addition of photon counting allows for univer- 
sal quantum computation. 

Indeed, CV clusters of up to four qumodes have already 
been experimentally realized [14, 15, 16J. In addition, 
recent results show that the network of linear optical el- 
ements may be eliminated entirely in favor of frequency- 
encoded qumodes and a single OPO [TH] EJ E2] • Such a 
method would be able to create a CV cluster state in just 
one step and in a single beam of light. Some such propos- 
als also have significant scaling potential 112) . What 
follows, however, will focus on the method described in 
Item 1 and discussed in detail in Ref. [9]. Item 2 sug- 
gests that once such clusters are available, they can be 
immediately tested by implementing protocols involving 
information distribution and other Gaussian operations. 
For example, this result immediately offers an experimen- 
tally viable method to use offline squeezed resources to 
perform squeezing operations online; such an online CV 
gate operation using offline CV resource states can then 
be not only "universally" applied to arbitrary optical sig- 
nal states [531 [2U [35], but would also no longer require 
adjustment of the offline resources to achieve different 
squeezing gates, as the CV cluster states provide a uni- 
versal resource for Gaussian computation together with 
homodyne detectors [39 . Finally, while accurate photon 
counting remains experimentally challenging, Item 3 im- 
plies that universal quantum computation is nevertheless 
possible. 



A. Cluster State Generation 

The naive canonical method to generate a given CV 
cluster state would be to apply the theoretical defini- 
tion directly, i.e., apply Cz interactions to a collection 
of squeezed states. While this method is conceptually 
simple, it is not very practical. The Cz operation does 
not conserve photon number, and requires the use of two 
single-mode online squeezers [40 . 




Linear Operator 




FIG. 4: A CV graph state can be generated by a Gaus- 
sian unitary acting on a collection of vacuum states [9]. The 
Heisenberg action of this Gaussian is given by a symplectic 
linear operator L acting on the phase space of quadrature 
operators (gi, q%, . . . ,pi,pi, ■ ■ ■)■ This linear action can al- 
ways be decomposed into the passing of offline squeezed states 
through a network of passive linear optical elements |40) . 



In a more practical approach, in Ref. [9 , it was shown 
that online squeezers are not needed at all. Any desired 
CV graph state of accuracy s is a pure multi-mode Gaus- 
sian state, and hence the only necessary online compo- 
nents are passive linear optics [HI HO] ■ To make this pre- 
cise, consider the generation of a graph state |</>(s)) corre- 
sponding to some graph G. Recall that \4>{s)) is defined 
by application of an appropriate sequence of Cz gates to 
a collection of squeezed states S(s) |0). 

The sequence of Gaussian transformations that take a 
collection of vacuum states to \<fi(s)) is represented suc- 
cinctly in the Heisenberg picture. Let v denote the vector 
of quadrature operators: v = {jq\, fe, ■ ■ ■ ,Px,p2, ■ ■ •)• I n 
the Heisenberg picture, the quadratures are transformed 
according to: 



v^CS(s)v = .M(s)v, 



(33) 



where S(s) represents the squeezing of each vacuum 
mode to form the state S(s) |0), and C represents appli- 
cation of Cz operations in accordance with the desired 
graph. Mathematically, these operations can be defined 
by their action on the quadrature operators: 



S(s)qi = scji, 
S(s)pi = pi/s, 



Cq t = qi 
Cpi = Pi 



(34) 



Concatenation of the two operations gives an explicit 
form for Ai(s). We refer to A4(s) as the generation ma- 
trix for \(f)(s)), which defines the Gaussian operation that 
generates |<^(s)) from the vacuum. The singular value 
decomposition for this matrix then provides an explicit 
recipe for how it may be generated using only linear op- 
tics and offline squeezing [3D] (See Fig. 4]). We refer to 
this method as the decompositional method 4 . 



4 In Ref. [9], the term "canonical" cluster states was reserved for 
those states that are obtained by directly applying a network 
of Cz gates onto momentum-squeezed states. These canonical 
states can then also be created, after Bloch-Messiah decomposi- 
tion, with offline squeezing and linear optics. In this sense, the 
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1. A simple example 

We illustrate the basic principles of this process by 
considering the explicit generation of the two-mode graph 
state, with generation matrix 



M(s) 



( s 
s 



s 
s 



(35) 



The four singular values of this matrix are given by Ai 
A2 = A+, A3 = A4 = A_ where 



\/\ + 2s 4 ±y/l + As 8 
V2s 



(36) 



A± specifies the amount of offline squeezing required to 
generate the two-mode graph state to an accuracy s (no- 
tice that A_ = A+ 1 ). That is, two squeezed states of 
magnitude A + (i.e., S(A+) |0)), together with passive lin- 
ear optical elements, may be used to generate this simple 
cluster state. Since we have transferred the online squeez- 
ing of the Cz operation into extra squeezing during the 
preparation process, A+ is generally greater than s. To 
generate the two-mode cluster to the same accuracy, our 
initial resources must be squeezed to a greater extent. In 
the case of the two-qumodc cluster, 



A. 



V2s 



s > 1. 



(37) 



Thus in the usual case where high accuracy is desired, we 
need to begin with states with a factor \/2 more squeezing 
to achieve a graph state of the same accuracy. This factor 
is known as the squeezing overhead. 

Since online squeezing generally represents a much 
greater experimental challenge than its offline equivalent, 
the decompositional method has a clear advantage over 
the canonical approach [9]. However, it is also fair to ask 



whether this method is superior in all situations. To test 
this, we consider the limiting case where offline squeez- 
ing is assumed to be as costly as online squeezing. One 
answer to this (in the affirmative) was already given in 
Ref. [5], wherein it is shown that extra local squeezing 
is required to obtain a canonically generated CV clus- 
ter state from an iV-mode Gaussian state in standard 
form [41] . We will revisit this result from another angle. 
In this case, the measure of the resource requirements is 
the total amount of squeezing required — whether online 
or offline measured additivcly in units of dB. 

We consider the toy case of the two-mode cluster here 
and follow with the general case in the following section. 
To generate this state up to accuracy s, the canonical 
method has two actions where squeezers are required: 

(a) Squeezing two vacuum states to S(s) |0), which re- 

quires two squeezers 101og(s 2 ) dB each. 

(b) Application of a single Cz gate (QND interaction), 

which requires two online squeezers of 4.18 dB 
each gD]. 

In contrast, the decompositional method requires squeez- 
ing in the following two steps: 

(a) Squeezing two vacuum states to S(s) |0), as before. 

(b) Squeezing these two states further by 101og(2) ~ 3 

dB each in order to account for the required squeez- 
ing overhead of v2- 

Thus the decompositional method saves (2 modes) x 
(4.18 dB/mode — 3 dB/mode) = 2.36 dB of squeezing 
for all values of s. In the next section, we show that 
the superiority of this method extends to general cluster 
states. 



2. Resource Requirements for General Graph States 



decompositional scheme is equivalent to the "canonical" scheme 
based on online Cz gates. In addition to those schemes resulting 
in the canonical cluster states, in Ref. [3], an alternative protocol 
was derived (independent of the Bloch-Messiah reduction), where 
offline squeezed states are sent through passive linear optics un- 
der the constraint that the outgoing multi-mode state satisfies 
the quadrature milliner conditions in the limit of infinite squeez- 
ing. In this case, for finite squeezing, there exist output states 
different from the canonical cluster states. This larger family 
of cluster states was referred to as cluster-type states, including 
many "non-canonical" cluster states, nonetheless satisfying the 
milliner conditions in the limit of infinite squeezing. With regard 
to experiments, an important feature of these generalized states 
is that the anti-squeezing components are suppressed by con- 
struction |I5I . In the present paper, in order to make a compar- 
ison between the scheme based on Bloch-Messiah reduction and 
that using direct Cz gates, the former is here referred to as the 
"decompositional method", while the latter shall be named the 
"canonical method" . Later, in order to make this comparison, 
no distinction will be made between offline and online squeezing. 



The above example motivates a general question: how 
much offline squeezing is required to create a graph state 
to accuracy s? In the case where Cz gates are directly 
applied, we would need n squeezed states of magnitude s. 
To generate such a state entirely by offline squeezing 
up to equal accuracy, the initial nodes would necessar- 
ily need to feature greater squeezing. The practicality 
of the decompositional method hinges on the size of this 
overhead. 

In this section, we show that there exist classes of uni- 
versal cluster states whose squeezing overhead per mode 
does not increase with the size of the cluster. In addi- 
tion, as in the case of two-mode clusters, the decompo- 
sitional method remains superior even if squeezing arbi- 
trary states (online squeezing) was as easy as squeezing 
the vacuum (offline squeezing). This is facilitated by a 
succinct method that computes the offline squeezing re- 
quired to generate any given graph state. 
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Theorem 2 We may generate an n-mode graph state 
with graph G to an accuracy s by passing n squeezed 
vacuum states — S'(sj) |0), with i — l,.,.,n — through a 
network of linear optical gates. Let A be the adjacency 
matrix of G I28f . In the limit of large squeezing fs> I), 
the level of squeezing required for each mode, Si, depends 
linearly on s, such that 



1 



(38) 



where ki are the singular values of A.. Thus, the squeezing 
overhead is bounded above by i/l + kf . 

Proof: Consider a particular graph state with generation 
matrix Ai(s). When the 1/s terms can be neglected, the 
n largest singular values of A4(s) then coincide with the 
singular values of the block-column matrix 



A 



(39) 



where I„ denotes the n x n identity matrix. The result 
follows. ■ 

Thus, the squeezing overhead of a given graph state is 
determined completely by the structure of its underlying 
graph. This theorem allows us to compute the exact 
resources required to generate any graph state. In certain 
situations, it is sufficient to know an upper bound on the 
amount of squeezed resources required. 

Theorem 3 To generate an n-mode graph state with 
graph G to an accuracy s (s 3> I), the maximal amount 
of squeezing required for any individual resource mode is 
bounded above by (Ks), where 



(40) 



K = a/1 + maxdeg z (G), 



and where maxdeg(G) denotes the maximum degree of G, 
i.e., the maximum number of edges connected to a single 
vertex. 

Proof: Consider first the case where G is an m-regular 
graph (i.e., each of its vertices has degree m). Let A be 
its adjacency matrix. Each row and column of A then 
sums to K. Therefore A may be written as the sum 
of m permutation matrices, P\, P2, ■ ■ . , P m - Noting that 
the largest singular value of a given matrix is its spectral 
norm (denoted ||-||), we have 



k, < HAH < 



3=1 



IP II <m, 



(41) 



in agreement with Eq. ( 40 1 . To generalize this result to 
an arbitrary graph state Grwith adjacency matrix A' and 
maximum degree m, observe that any graph of maximum 
degree m may be obtained by removing edges from an 
m-regular graph G. The spectral norm of an adjacency 
matrix strictly decreases with the removal of an edge, so 
we have ||A'|| < ||A|| < m. ■ 



This theorem immediately implies that the squeezing 
overhead for universal cluster states of any fixed accu- 
racy is bounded. Such states have maximal degree of 4, 
and hence feature a squeezing overhead of VT7. Thus, to 
guarantee the generation of a universal cluster to accu- 
racy s, one would need to (a) generate a lattice of optical 
modes, each of which is squeezed up to 10 log (s 2 ) dB 
and (b) proceed to squeeze each mode by a further 12.31 
dB. Meanwhile, quantum wires would require a maxi- 
mum overhead of a/5, and hence an additional 6.99 dB 
of squeezing. 

To see that the decompositional method is superior to 
the canonical method, recall that each Cz gate requires 
two 4.18 dB squeezers. In the case of a universal square 
lattice, the ratio of edges to vertices is 2:1. Thus, while 
the decompositional approach requires at most an addi- 
tional 12.31 dB per vertex, the two Cz gates applied per 
vertex would cost 16.72 dB (since each Cz gate requires 
two online squeezers of magnitude 4.18 dB). For a square 
lattice of size N x TV, the decompositional method would 
save approximately 4.417V 2 dB of squeezing. Thus, not 
only is it more practical to perform universal quantum 
computation through squeezed offline resources, but it 
also turns out to be more efficient, i.e., cheaper in terms 
of squeezing resources required. 

A typical setup would involve the generation of n 
squeezed optical modes. These are then passed through a 
network of linear optical gates, of which the resulting en- 
tangled beams formally encode the desired cluster. Since 
there exist cluster states that are universal, the setup of 
this generation process does not need to be altered for 
different algorithms and hence may potentially be mass 
produced. The resulting beams can then be measured to 
perform the desired quantum computation. 

As a final remark in this section, we come back to the 
question whether redundant, multiple-rail encoding may 
suppress the accumulation of finite-squeezing errors in a 
cluster computation [5] . According to Eq. ( 40 1 , the offline 



squeezing per mode for generating a multiple-rail graph 
of accuracy s with m rails [5] is bounded above by (Ks) 
with K = \/l + fn 2 ■ Therefore the initial squeezing vari- 
ances may be as small as l/(2s 2 ) x 1/(1+to 2 ), so roughly 
l/(2s 2 m 2 ) for large m. This lower bound converges to 
zero faster than the actual reduction of the excess noise 
in the cluster computation which scales as l/(2s 2 m). If 
one has access to squeezing resources with variances of 
l/(2s 2 m 2 ), one may better use them directly without 
multiple-rail encoding [9]. However, note that also here 
no complete proof for the effective failure of a decomposi- 
tional multiple-rail scheme is given, as the above analysis 
only relies upon bounds. 



B. Optical Cluster-State Computation 

The measurement of optical modes completes the op- 
tical implementation of CV cluster-state computation. 



Recall from Section IV B that any multi-qumode oper- 
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ation may be implemented by measurements that gener- 
ate (a) the shearing transformation e ls ^ I" 1 and (b) the 
cubic phase gate e ls ^ ^ (while a nontrivial two-mode 
gate only requires measuring the p quadratures on a two- 
dimensional cluster state). 

The implementation of (a) is reasonably straightfor- 
ward. The required measurement basis, p + sq, is a ro- 
tated quadrature basis. We can write it in the standard 
form r[sm(8)q + cos(8)p] — rpg, where 



tan(0) 



(42) 



Thus, the optical implementation involves measurement 
in the rotated quadrature basis pg, followed by a rescal- 
ing of the result by a factor of r. Therefore, all multi- 
qumode Gaussian operations may be achieved via sim- 
ple quadrature measurements on a sufficiently connected 
graph state. While such operations are insufficient for 
universal quantum computation, they allow for general 
graph transformations, and thus many optical experi- 
ments that test the foundations of cluster-state quan- 
tum computation can be performed using just quadrature 
measurements — i.e., homodyne detection. 

The physical implementation of the cubic phase gate 
is more challenging. Since the required Hamiltonian to 
be implemented is no longer quadratic in the quadra- 
ture operators, a nonlinear optical element is required [1J. 
Two separate strategies may be employed, involving ei- 
ther embedding the nonlinear resource within the cluster 
(making it a non-Gaussian state) or using photon count- 
ing on the modes of an existing (Gaussian) graph state. 



1. Quantum Computation by Nonlinear Resources 

In the standard model of cluster-state computation, 
all qumodes are initialized in the state |0) prior to the 
entangling operation. The quantum computation is en- 
tirely dictated by our choice of measurement basis D^pD. 
However, it is also possible to encode the computation- 
gates D within our initial resource. To see this, note that 
since D commutes with Cz, the circuit 



|0> P ■ 

is operationally equivalent to 



DipD )= 



771l 



777,2 



D 



P )= 7771 



-jjT)= m 2 



(43) 



Therefore, instead of measuring in the D^pD basis, we 
may have used instead the resource state D |Q) for cre- 
ating the initial cluster. Since the cubic phase gate 



D = e ls * 3 / 3 allows for universal quantum computation, 
this observation suggests that the cubic phase state, 
e ls <7 /3 \0) p , will have the same effect. One method to 
optically generate such states is given in Ref. 

Should these states be available, they may be attached 
at set locations within a universal cluster state in place of 
the usual |0) . The resulting non-Gaussian cluster would 
be an improved resource for universal quantum compu- 
tation. Any CV unitary may be implemented employ- 
ing such clusters, even when one is limited to quadrature 
measurements only. Of course, generating non-Gaussian 
quantum states remains an experimental challenge, and 
the cubic phase state is no exception. 



2. Cluster-State Implementation of the Cubic Phase Gate 

The previous observation suggests that if there exists 
a method of generating the cubic phase state by viable 
single-qumode measurements on a standard cluster state, 
then cubic phase gates may be applied to arbitrary in- 
puts. 

One possible approach [4 involves photon counting on 
one arm of a displaced two- mode squeezed state. This 
procedure may be summarized by the following quantum 
circuit: 



5(010) 
S(s) |0) 









Z{r) 




B 





77 )= n 



e i-i(n)q 3 



(44) 



This circuit entangles two highly squeezed states, S(s) |0) 
and Sis' 1 ) |0), with s > 1, via a standard beam- 
splitter interaction B. A large momentum displace- 
ment Z(r), with r > s, is applied to the resulting two- 
mode squeezed state. We then make a photon count- 
ing measurement 77 on the displaced mode, which ap- 
proximately collapses the unmeasured qumode into a 
cubic phase state e* 7( -™^ |0) , dependent of the mea- 
surement result 77 through 7(77) = (6-y/277 + This 
procedure is essentially a measurement-based quantum 
computation — it involves the application of suitable mea- 
surements on an entangled resource. Thus, we may recast 
it into the form of a standard cluster-state computation. 
The two-mode squeezed state generated coincides with a 
two-qumode cluster state modulo a local Fourier trans- 
form on one of the nodes. These observations allow us 
to construct a circuit that is consistent with the cluster- 
state formalism, while also functionally equivalent to Cir- 



cuit (441, i.e., 



S(s) |0) 
S(s) |0> 





X(r) 






1 


1 





(45) 



In this circuit, the initial entangled resource is a standard 
two-qumode cluster arranged in a linear configuration. 
The quadratures of this state are rotated with respect to 
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the preparation using a simple beamsplitter, so a position 
displacement is applied to the first qumode, followed by 



a photon counting measurement. Just as in Circuit (44 1 



the second qumode is then collapsed to an approximate 
measurement-dependent cubic phase state. 

Should the above circuit be attached to the second 



qumode of Circuit ( 43 1 , we may apply the cubic phase 
gate e* 7 '™^ to an arbitrary input state. A cubic phase 
gate e laq , for any a, may be decomposed into a combina- 
tion of e i7 (™)9 3 and two squeezers that depend on both n 
and a [I]: 



S\t(n))e l ^ n) i 3 S(t(n)) = e la *\ t(n) = [a/j(n)] 



1/3 



The addition of these squeezers to Circuit ( 43 1 gives 
a measurement-based scheme to implement e laq for 
any desired a (modulo measurement dependent shifts in 
phase space): 



S(t(n)) 



JL 



mi 
m 2 



(46) 



SHt(n)) - e^ 3 |0) p 



Since the squeezing operation S(t(ri)) is Gaussian, it can 
always be implemented by a suitable sequence of quadra- 
ture measurements on a linear cluster state. Combining 
this circuit with Circuit ( |45| leads to a cluster-state im- 
plementation of the cubic phase gate (see Fig. [5j. 

Observe that since the squeezing strength t(n) is de- 
pendent on the outcome n of the photon counting mea- 
surement (in addition to the fixed parameter a), imple- 
mentation of the squeezers must be done after the photon 
counting. As with non-Clifford group operations on qubit 
clusters, adaptive measurements are involved, and hence 
the order of the measurements now matters. 



VIII. DISCUSSION AND CONCLUSION 

Quantum information processing through the use of 
CV cluster states is a recent development in the field of 
quantum computation. While its basic principles have 
already been introduced in Ref. |7J, here we fleshed out 
the details of the protocol in a way that we hope will fa- 
cilitate further understanding of CV cluster states, along 
with their potential optical implementation. 

The Schrddinger representation of CV graph states 
used in Ref. [7] becomes unwieldy as we explore graph 
states of nontrivial size. In such cases, Heisenberg mil- 
liners and Wigner functions are potentially useful tools. 
We showed how these representations may be utilized 
to derive rules of thumb on how graph states transform 
under measurements, as well as the existence of a CV 
cluster state that can be used as a resource for arbitrary 
CV operations. 
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FIG. 5: A proposed cluster-state implementation of the cu- 
bic phase gate e laq applied to an arbitrary input state \<f>). 
The displaced number-state measurement on node C imple- 



ments Circuit ( 45 1 , resulting in the generation of a cubic phase 



state at D. This nonlinear resource state, together with sub- 
clusters that generate the n-dependent squeezing corrections 
(B and E), allows us to apply a cubic phase gate of any spec- 
ified strength. 



The optical implementation of CV cluster states has 
also been further explored. When the decompositional 
method was initially introduced in Ref. [5], one of the 
primary concerns was that the price we pay for avoid- 
ing the need to perform online squeezing was an exces- 
sively large overhead in the extra offline squeezing re- 
quired. Our results alleviate this concern. We proved 
that the squeezing overhead per mode, when only offline 
squeezing is used for a given accuracy of the CV cluster, 
does not increase with the size of a universal CV clus- 
ter state. The upper bounds derived on the necessary 
amount of offline squeezing indicate that the decompo- 
sitional approach has significant advantages over the di- 
rect approach through QND interactions [7J, even in cases 
where online squeezing is no more costly than its offline 
counterpart. While universal quantum computation us- 
ing CV cluster states may be no less challenging than its 
qubit counterpart, the generation of CV clusters- either 
using the decompositional method |Hj or through one-step 
generation using a single OPO |llj — is potentially more 
viable than in corresponding optical qubit schemes [8]. 

To perform universal quantum computation, we 
adapted the experimental generation of the cubic phase 
gate as given in Ref. [I] to the cluster-state formalism. In 
addition, by extending the CV cluster-state framework 
to include the use of non-Gaussian resource states, we 
showed that possession of a suitable non-Gaussian state is 
sufficient for universal quantum computation within the 
cluster-state framework, even when only Gaussian opera- 
tions, i.e., homodyne detections are employed during the 
cluster computation. This leads to promising possibilities 
for universal CV cluster-state computation. One could, 
for example, envision that difficult nonlinear measure- 
ments are used to generate non-Gaussian resource states 
offline, which may then be distributed to consumers who 
are limited only to simpler measurements — i.e., quadra- 
ture homodyne detections. The consumers can neverthe- 
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less use the non-Gaussian states as resources for universal 
quantum computation. 

In theory, all the ingredients for universal CV cluster- 
state computation have been developed. In particular, 
since the necessary squeezing resources for creating a 
cluster of given accuracy have been shown to be inde- 
pendent of the size of a universal cluster state, scalability 
would only depend on the ability of suppressing the accu- 
mulation of errors at every measurement step during the 
cluster computation. At least for homodyne detections 
with near-unit efficiency, these errors mainly originate 
from the finite squeezing and they grow linearly with 
the length of the cluster. This could be compensated 
by increasing the accuracy of the cluster, hence making 
the squeezing per mode again dependent of the size of 
the cluster and the computation. Alternatively, some 
form of error correction may achieve full scalability in a 
strict sense, similar to fault tolerant schemes for qubit 
quantum computation. To find an efficient method for 
CV quantum error correction, combined with a practical 
scheme for incorporating a non-Gaussian element into 
the Gaussian CV cluster-state framework, remains the 
main challenge to scalability and universality of the CV 
approach. Nonetheless, our results here are an impor- 
tant step towards small-scale proof-of-principle demon- 
strations of CV cluster computation. 
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APPENDIX A: THE NULLIFIER FORMALISM 
FOR QUADRATURE MEASUREMENTS 

Graph state transformations through quadrature mea- 
surements have an efficient nullifier description. Consider 
a measurement pi on a graph state \<j>) nullified by a vec- 



tor space with basis elements {Hi, H 2 , . ■ . , H n }. There 
are only two distinct cases: 

1. pi commutes with all basis elements Hj. Then, \<j>) 
must be an eigenstate of pi with some eigenvalue 
m l . 

2. There exists exactly one basis element that does 
not commute with pi ; call this clement Hi . 

This is because if there were to exist two or more basis 
elements that do not commute with pi, say Hi and Hj, 
then there exists a constant k such that [pi, Hi + kHj] = 
0. Thus, it is always possible to construct a basis such 
that only one element does not commute with pi, with 
\pi,Hi] = -i. 

In case (1), H^pi \4>) — piHk \<f>) — for each nullifier, 
and thus pi\4>) — mi \<p) for some m^. Hence {pi — mi) \4>) 
must be a nullifier of \<p) for some m,. The measurement 
of pi yields rrii, and the state remain undisturbed. 

In case (2), we write Hi as qi + J2j £jPj + djfa + 
Co for some constants Cj and dj. Let the measurement 
result be m^. Then, the transformed nullifier algebra is 
obtained by replacing Hi with pi — mi. Since the i th 
mode is now disentangled from the cluster and no longer 
interesting, we may discard it by choosing a basis such 
that all but one of the elements (namely pi — m t ) acts 
as the identity on the i th mode. Measurement in the 
computational basis q can be analyzed analogously. More 
general quadrature measurements of the form p s ^2 j 2 — 
p + sq can be treated in this formalism by application 
of the unitary e ls ^ I" 1 , followed by a standard momentum 
measurement. 

We demonstrate this formalism on the case where p- 
measurements are made on the first two qumodes (with 
measurement results m\ and m-i) of a linear three- 
qumode cluster, defined by the three milliners 

{H X ,H 2 ,H Z } = {pi - fa,p 2 ~q\- fa,p3 - fa}- (Al) 

pi commutes with all milliners except H 2 , so we replace 
H 2 with pi — mi, giving the milliners {pi — fa, pi — m,p^ — 
fa}, which defines the same state as {fa — m\,pz — m\} 
after discarding the measured mode. 

Repeating this procedure for the measurement of p 2 
results in the nullifier of the output state being p^ — m\. 
Thus the remaining unmeasured node is in the state 



\mi) p , in agreement with Eq. (15 1 
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